Diffraction of light by a planar aperture in a metallic screen 
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o 

. We present a complete derivation of the formula of Smythe [Phys. Rev. 72, 1066 (1947)] giving 

' the electromagnetic field diffracted by an aperture created in a perfectly conducting plane surface. 



The reasoning, valid for any excitating field and any hole shape, makes use only of the free scalar 
Green function for the Helmoltz equation without any reference to a Green dyadic formalism. We 



g . compare our proof with the one previously given by Jackson and connect our reasoning to the general 

Huygens Fresnel theorem. 

P ■ PACS numbers: 

o 

p I I. INTRODUCTION 



Diffraction of electromagnetic waves by an aperture in a perfect metallic plane is not only a mathematical problem 
of fundamental interest but is connected to many applications in the microwave domain (for example, in waveguides 

The 



and in cavity resonators [1| ) as well as in the optical regime where it is involved in many optical arrangements 
fNl fundamental importance of this phenomenon in near-field optics has been pointed out as early as in 1928 by Synge 

o . . . 

in his prophetic paper and is currently involved in modern near-field scanning optical microscopy (NSOM) 



In the domain of applicability of NSOM where distances and dimensions are smaller than or close to the wavelength 
of light, we need to know the exact structure of the electromagnetic field, and we cannot in general consider the 



usual approximations involved in KirchhofF's theory for a scalar wave 
approaches is the one given by Bethe 



5|-|7[. In this context, one of the most cited 



8| in 1944 and corrected by Bouwkamp P, [l^. It gives the electromagnetic 
field diffracted by a small circular aperture in a perfect metallic plane in the limit where the optical wavelength is 
much larger than the aperture. Less known is the more general formula of Smythe [ill , [l^ which expresses in a 
formal way the Huygens Fresnel principle for any kind of aperture in a metallic screen. Even if this formula is not 
an explicit solution for the general diffraction problem, it constitutes an integral equation which can be used in a 
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self consistent way in perturbative or numerical calculations of the diffracted field [13|, [IJ] . Further efforts have been 



made by Smythe 
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HQ. 



12| himself in order to justify his formula by means of some arrangements of current sheets 
fitting the aperture. This method essentially consists of transforming the problem of diffraction by a hole into a 
physically different one in order to guess the correct integral equation for the original problem. However, if this 
physical reasoning proves the consistency of the proposed solution with Maxwell equations and boundary conditions 
for the field, it is not directly connected to the rigorous electromagnetic formulation of the Huygens Fresnel principle 
obtained by Stratton and Chu [l^. Such a connection is expected naturally because these two formulations of 
diffraction must be equivalent here. 

Jackson |lgf|, in the first edition of his textbook on electrodynamics, developed a complete proof of the Smythe 
formula starting from the Stratton and Chu formula [Eq. ([3|) of the present paper]. Nevertheless, like in the original 
paper of Smythe, Jackson transforms the problem into a physically different one in order to guess the correct result. 
The result is then subjected to the same remarks as above for Smythe's approach. Other justifications of Smythe 
results are based on the use of the Babinet theorem or of the Green dyadic method. The latter, which uses a tensorial 
Green function instead of a scalar one like in Kirchhoff's or Stratton and Chu's theories, gives us the most direct 
justification for Smythe approach in terms of the Huygens Fresnel principle. However, this proof is for the moment 
not directly connected to the Stratton and Chu approach. It is the aim of this paper to establish such a link. 
The paper is organized as follows. We give in Sec. H a description of the general theory of diffraction of electromagnetic 
waves by an aperture in a screen. In Sec. HI, we exploit precedent works by Jackson [igIIitIi and Levine and Schwinger 
[isl to justify directly and rigorously the Smythe formula using the Stratton Chu theorem without relying on any 
ingenious physical "trick" . Sec. IV deals with a vectorial justification of Smythe's approach. The consistency between 
the various theoretical treatments of diffraction by an aperture in a metallic screen is stressed in Sec. V which also 
compares our treatment with that obtained within the Green dyadic formalism [igl . [2^ . Our conclusions appear in 
Section VI. 



II. THE DIFFRACTION PROBLEM IN ELECTROMAGNETISM 

The first coherent theory of diffraction was elaborated by Kirchhoff (1882) on the basis of the Huygens Fresnel 
principle [2, 21 1. The method of integral equations allows one to write a solution ?/;(r^e^'"* of the Helmholtz 
propagation equation [V^ + k^]ijj (f) = {k = lo/c) using the "free" scalar Green function G |^r, j''^ — e^^^ /AttR which 
is a solution of the equation [V^ + k^]G ^r, r'^ = —5'^ if ^ ''') ' 

If, as schematized in Fig. 1, we consider now an aperture 5S made in a two-dimensional infinite screen S and 
illuminated by incident radiation, we can express the field tp existing at each observation point located behind the 
screen (i. e. , for z > 0) by the Kirchhoff formula 

1/, (f) = / [Vj (r'j n' ■ V'G [r, r') - G (r, n' ■ V'^j (r'jjdS", (1) 



where the normal unit vector n' is oriented into the diffraction half-space. 

In a problem of diffraction, we usually impose the additional first KirchhofF "shadow" approximation ^r'^ = 
dn''4' ij"'^ ~ ^ which is valid on the unilluminated side of the screen. This permits one to restrict the integral in ([1]) 
to the region of the aperture only, which is very useful in some approximations or iterative resolutions. Nevertheless, 
this intuitive hypothesis has some fundamental inconsistencies because, following a theorem due to Poincare 21|, a 
field satisfying the shadow approximation on a finite domain must vanish everywhere. 

A classic solution proposed by Rayleigh |22'] and Sommerfeld 23] to circumvent this difficulty consists in replacing 
the free Green function by the Dirichlet Gd or the Neumann Gn Green functions [l^ satisfying dn'GN (j^,"!"'^ = 
and Gd (^t^, ^'^ ~ for all points r' on S. We can then rigourously reduce the integral to the region of the aperture 
depending on the nature of the boundary problem. For example, if we impose ■0 = on the screen, we can then write 

(r) = / (?) dn'Go (r, P) dS'. (2) 

J Aperture 

In principle, it could be possible to generalize the preceding methods to the different Cartesian components ipa of 
the electromagnetic field using equations of the form i/ja = Jgii/Jadn'G — Gdn'ipa\dS' . Nevertheless, as pointed out 
by Stratton, Chu [15] and others 
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-26|, the Maxwell equations couple the field components between them and 
the consistency of these relations must be controlled if we use an integral equation like Eq. ([T]) either in an exact or 
approximative treatment of diffraction. In addition, because the boundary conditions imposed by Maxwell's equations 
connect the tangential and the normal components of the field on the screen surface, it is not at all trivial to reduce 
the integral to the region of the aperture directly using Eq. 

Due to the uniqueness theorem, such possible reduction of the integral appearing in the Huygens Fresnel principle is 
expected in the case of a perfectly conducting metallic screen. Indeed, following this uniqueness theorem, the field in 
the diffracted space must depend only on the tangential electric field on the screen and aperture surface. Because the 
tangential electric field vanishes on the screen, the integral must depend only on the tangential field at the opening. 
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Numerous authors, especially Stratton and Chu [l^ as well as Schelkunoff 
equation satisfying Maxwell's equations automatically. We can effectively write 

1^ {x)^ ( [ik {n' x1^)g+ (n' xl^) x V'G - 
J s 



281 ] , have discussed a vectorial integral 



(n' ■ ~t) V^G]dS", (3) 



hereafter referred to as the Stratton Chu equation. A similar expression holds for the magnetic field by means of the 
substitution ^ it and it ^ -it. 

It is important to note that Eq. ^ is over-determined although it depends explicitly on the tangential and normal 
components of the electromagnetic field defined on 5*. Indeed, due to the equivalence principle of Love and Schelkunoff 
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29| and to the uniqueness theorem, we expect that the "most adapted" integral equations depend only on 
n' X 1^ or n' X 1^ on S. In addition, unlike in the scalar case, we cannot directly reduce the surface integral to the 
region of the aperture just by choosing an adapted Dirichlet or Neumann Green function. It seems then necessary to 



apply once again the shadow approximation of Kirchhoff in order to simphfy the integration despite the inconsistency 
of the method. As in the Poincare theorem, some problems appear here because we need to add a nonphysical 
contour integral associated with a magnetic line charge in Eq. ([3]) (or to an electric line charge in the equivalent 
formula for B) in order to satisfy Maxwell's equations and to compensate for the arbitrary change imposed to the 



integration domain [3^. Furthermore, in this Kirchhoff Kottler [2(| theory, the introduction of contour integrals 
induces a logarithmic divergence of the energy at the rim of the aperture, a fact which is forbidden in a diffraction 
problem. 

The particular case of the diffraction by an aperture in a planar screen constitutes an exception in the sense that 



a rigorous integral equation had been anticipated by Schelkunoff '27] and Bethe 



8| for a subwavelength circular 



aperture and generalized by Smythe [ill [l^ for any kind of aperture. The integral equation is 

For some applications, it is important to note that in the short wavelength limit (A ^ aperture typical radius) for 
which the electromagnetic field in the aperture can be identified with the incident plane wave (first 
Kirchhoff approximation), the formula of Stratton Chu limited to the aperture domain and the exact solution of 
Smythe give approximately the same result. Indeed, within the Fraunhofer approximation, Eq. Q reads 



i^^^fx/ fi4^e---^)d5', (5) 



Aperture 



2tt 



whereas Eq. ([3|) reduces to 



^ 2 ./Aperture ^ ^TT 



Both equations are identical in the practical limit of small diffraction angles, i. e. , close to the normal axis z going 
through the aperture. Equation ([5]) is correct for a subwalength aperture only because we cannot identify the field in 
the aperture with the incident one. We can see that the asymptotic diffracted field for z > is equivalent to the one 
produced by an effective magnetic dipole 

iS..^/ (7) 



Aperture 



2'Kik 



and by an effective electric dipole 



^cff = ^ / (x' ■ dS'. (8) 

4"" v/ Aperture ^ ^ 

These formula are fundamental in the context of NSOM because they give us the Bethe Bouwkamp [^-Q Q 

dipoles 

which, in the particular case of a circular aperture of radius a, are 

^c^fi^f M=-'-ft\r (9) 

OTT OTT I 

and ^\"^ are, respectively, the locally uniform normal electric field and tangential magnetic field existing in the 
aperture zone in the absence of the opening (in z = 0^). 
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III. GREEN DYADIC JUSTIFICATION OF THE SMYTHE FORMULA 

The so-called Smythe formula Eq. (j?]) is generally obtained on the basis of different principles such as the Babinet 
principle or the equivalence theorem (see Schelkunoff Bouwkamp [30], Jackson [l7|). In particular, the 

equivalence theorem shows that the solution of Smythe for z > is identical to the one obtained by considering a 
virtual surface magnetic-current density given by = -cz X t/ (27r) . All these derivations are self consistent if we 
consider the very fact that the guessed results fulfill Maxwell equations. Then, the uniqueness theorem ensures that 
the result is the only one possible. Nevertheless, as already noted, the calculation is not direct and not necessarily 
connected to the Stratton and Chu formalism. A classical calculation due to Schwinger and Levine [l9|, |20[ shows. 



however, that it is possible to rigourously and directly obtain this equation using the tensorial, or dyadic. Green 
function formalism. 

. n « . . . 

Such an electric dyadic Green function [31!] Ge, which is solution of the equation 

^ X (^x Ge (r,/)) = fc2 Ge (f.r') + 1 5^^ (t - r'^ (10) 

(with 5 ~ X)i S:iXi) satisfying the condition Ge= — (l/^^) ^ ^ ' '^^'^ used to write the integral equation 

'^{r)^ I [{n' xl^yV'xGe ~ik^ ■ (?i'x Ge)]dS" (11) 
^ s 

which is defined on the same surface as previously. By imposing the dyadic Dirichlet condition n' x Ge— on 5, we 
can obtain the relation 

'^{r^= [ [(n' xt)-V'x Ge]dS' (12) 

J Aperture 

which depends only on the tangential electric field at the aperture. This is in perfect agreement with the equivalence 
principle and the uniqueness theorem. 

n 

Following Ref. [Slf, the total Green function Ge for the plane can be deduced from the "free" dyadic 



^0 
Ge 



[with R—^{x — x')^ + {y ^ y')^ + ^ ^0^] by using the image method. We have 



Ge (.>') = -^^V] Gn {ry)+2H'^, (14) 
where Gd — ^e'^^^/i? — e**^^ /An is the scalar Dirichlet Green function for the plane screen, and R' = 



[x — x')^ + {y ~ y')^ + {z + z')^. Inserting this Green function into Eq. (IT2]) gives us directly Eq. ([4]). It is in- 
teresting to observe that with the Green dyadic method, we can recover the formula of Smythe by using a magnetic 
current distribution located in front of a metallic plane or, equivalently, by using a double layer of magnetic currents 
propagating in the same direction [l^ • 



In theory, both approaches based either on the scalar Green functions or on the dyadic Green functions are equiva- 
lent. In practice however, the difficulties related to the Stratton Chu formula Eq. ([3]) have imposed the Green dyadic 
method. An illustration of this statement is that the dyadic formalism has been extensively used in the context of 
the electromagnetic theory of NSOM 331-136 1. 



IV. VECTORIAL JUSTIFICATION OF THE SMYTHE FORMULA 

We propose now a justification of Eq. (U) based on the Stratton Chu formula Eq. (jS]). This derivation will directly 
reveal the equivalence of the scalar and dyadic approaches in the particular case of a planar screen with an aperture. 
Let the surface S of equation z — be an infinite, perfectly conducting metallic screen containing an aperture covering 
the surface SS. By the definition of diffraction, we can always separate the total electric (magnetic) field into 
an incident field i^* (i^* ) existing independently of the presence of the screen, and into a diffracted field 1^' (i^') 
produced by the surface charge and current densities p^, J' s located on the metal. 



We have B' — xA' and E' = + ikA' where potentials are expressed in a Lorentz gauge 

./Screen V C V / R J 



$'(r1= / dS'(p'JrY—-), 

./Screen ^ \ / H / 



(15) 



with R = 11?^ — 7"' II (we omit here the time dependent factor e "^*). Because these potentials are even functions of z 
we then have the following symmetries 

E'^, E'y, B'^ are even in z, 

B;, are odd in z. (16) 

These symmetries already used by Jackson [l^,[l3| imply in particular E'^ — By = B'^ = at the aperture. Therefore, 
the field is a discontinuous function through the metal. 

Let us now consider an observation point x located in the half space z > 0. We can apply the vectorial Green theorem 
on a closed integration surface made up of a half sphere "at infinity" and of the S'^ plane {z = 0^) as seen in 
Fig. 2 (A). This surface 5*+ can itself be decomposed into an aperture region 6S^ and into a screen region (5* — SS)^ . 
We have then 

— > /■ / ^ — >\ / ^ — >\ — > / ^ — >\ — > f / ^ — >\ — > 



E' (x) = [ik (n' X B') G + (n' x E'^ x V'G + {n' ■ E'^ V'G]dS" + J [(n' x E'^ x V'G]dS" 

+ I [ik (n' xB'^G+in'xE'^ x VG + (n' ■ E'^ V^G]dS", 

(17) 



where the unit vector n' lies on S'^ and is oriented in the positive z direction: n' — z. Similarly we can consider the 
surface of integration represented in Fig. 2 (B). We obtain an integration on the 5+ , surfaces and on {S — 6S)^ 
and (5* — 5S)^ surfaces. Such integration surfaces have already been used by Schwinger and Levine in the context of 
diffraction by a scalar wave [l^. Here, due to the symmetries given by Eq. (fT6l) . we deduce 

'e' [x)^2 [ik {n' xB'^G+ (n' ■ i^j V'G]dS' + j [ik {n' x 13'^ G + (n' x E'^ x V'G + (n' ■ IS'^ V^G]dS" 



+ J [ik (n' xI3'^G+(n'xE'^x V'G + (n' ■ 1^'^ V^G]dS" 

(18) 



with n' — z on the [S — 58)'^ surface. After identification of Eq. ([T7)) and Eq. ([T^ . we obtain 



'E'[x)^2j [{n' xl^'^ x^'G]AS' - j [ik i^n' xIb'^ G + i^n' xl^'^ x^' G + i^n' -E'^^' G]AS' 

[ik [n' x~B'^G+ {n' xl^'^ x V'G + [n' ■ V^G]dS". 

(19) 

In order to simplify this formula, it is important to note that the fields i?', B' located on are the reflected fields 
which could be produced by the complete metallic screen z = submitted to the same incident field in the 
absence of the aperture. 

Because this field compensates for the incident field for z > 0, we have in this half-space. As 

a consequence, the integral on in Eq. (|19[) can be written 

[x) + [ik (n' x^'^G+ [n' x x V'G + (ij' ■ 1^'^ V^G]dS", which is a direct application of the Green 
theorem for an observation point located on the closed surface composed of and 5*^. 

Injecting this last result into Eq. and after subtracting and adding 2 Jg+[(ji' x i^'j x V'GJdS", we finally obtain 
E^' = + ii't^) ^ijg^g 

{x) = 2 J [(n' X i^) X V^G]d5' - 1^' {x) dS' (20) 

and 



(f) = _ y [ik X G + (n' X 1^'-^ X V'G + (n' ■ i?' ) V^G]dS". 
+ J [ik (^n' xl^'^G- (r? X i^') x V'G + (n^' ■ i^') V^G]d5'. 



(21) 
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Because of Eq. ([16]), we also have 



and 



E^,y (x, y, z) = -El y {x, y, ~z) , 
K (a;, y, z) = -B\ (x, y, ~z) , 

Bl,y{x,y,z) = Bly{x,y,-z), 
K {x,y,z) = El {x,y, -z) 

(22) 



for z < 0. Using the fact that the integral on S*"*" can be written as an integral on S~: = — jg- {^', ^'}, 

and using Eq. (|22|) , the last two integrals in Eq. \2\\ can be transformed into [ik {n' x i^''^ G+ ^n' x i^''^ x 

V'G + ■ is' '"^ V'GJdS". Because the observation point is outside of the closed surface composed of and of S*^, 
i^'-^^ (!?") is zero. Regrouping all terms, the total electric field in the half plane z > is finally given by the Smythe 
formula: 



i^(^)=2/ [fn' x i^) X V^G]dS" 
Jss+ ^ ' 



ikR 

; (23) 

^ J Aperture / n / 

where we have applied Maxwell's boundary conditions that annihilate the tangential component of the total electric 
field on a perfect metal. An equivalent derivation in the z < half space gives 



i^(l^) = i^°(^) + 2^ [(n' X i^) X V^G]dS" 



= ^" (-i^) - X ( / (z X '-^dS') , (24) 

V ^Aperture ^ J H J 

where is now the total electric field existing in the z < domain for the problem without aperture. 



V. CONSISTENCY BETWEEN VARIOUS APPROACHES 

As written in the introduction, the proof given by Jackson [l^ of the Smythe equation is connected to the theory 
of vectorial diffraction Eq. In order to solve the problem, Jackson used a volume looking like a flat pancake 
limited by the two and surfaces, and he applied Eq. ([3]) to this boundary. Then, in agreement with Smythe, 
Jackson imagined a double current sheet such that the surface current on the two 5*+ and layers at any point of a 
given area fitting the aperture are equal and opposite. With such a distribution, it is possible to reduce the integral 
of Eq. to the one given by the formula of Smythe, Eq. (1^^ . Such a formula is then the correct one to describe the 
diffraction problem by an aperture in agreement with the uniqueness theorem. 

Our justification of the Smythe theorem is more direct because it uses only the Huygens Fresnel theorem without 



applying the intuitive trick of a virtual surface current distribution associated with a different physical situation 
(double layer of electric current, or layer of magnetic current confined to the aperture zone). Our result is in fact the 
direct generalization of a method used by the authors for a scalar wave ip. Using two different surface integrations, as 
the ones used in this paper, we are indeed able to prove directly the Rayleigh-Sommerfeld theorem given by Eq. (p)). 
This scalar reasoning, which is similar to the one presented before, is given in the appendix. It can be observed that 
the scalar result makes only use of the Green function in vacuum G in order to justify the result obtained with the 
Dirichlet one Gd- Similarly, our derivation of the Smythe formula uses the scalar Green function in order to justify 
the result obtained with the "Dirichlet" dyadic Green function. Then, the two reasonings presented in this paper 
for an electromagnetic and a scalar wave show the primacy of the Huygens Fresnel theorem given by Eq. ([1]) for the 
scalar wave and by Eq. ([3|) for the electromagnetic field, respectively. 

A few further remarks are here relevant: First, the mathematical results described here constitutes a justification of 
the physical "trick" introduced by Smythe and Jackson. However more work have must be done in order to see if the 
method based on scalar Green functions could be extended to other geometries. Second, the Smythe formula allows 
one to express the electromagnetic field radiated by the aperture (far-field) as a function of the near-field existing in 
the aperture plane. This method could thus be useful for calculating the field generated by a NSOM a pert ure if we 
know the optical near-field (computed, for example, by using numerical methods discussed in Refs. 
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VI. CONCLUSION 



In this paper, we have justified the vectorial formula of Smythe expressing the diffracted field produced by an opening 
created in a perfectly metallic screen. Our justification is based only on the Huygens principle for electromagnetic 
wave and on the specifical nature of boundary conditions for the Maxwell field. This proof differs from the ones 
presented in the literature because it does not use the concept of current sheets introduced by Smythe and Jackson. 
The demonstration uses only the scalar Green function in free space and does not consider Dirichlet or Neumann 
boundary conditions as involved in the Green dyadic method. 



Appendix A 

Let 5* (r) be a scalar wave solution of the Helmoltz equation for the problem of diffraction by an opening SS in 
a plane screen S. In order to define completely the problem, we must impose boundary conditions on the screen 
surface. Here, we choose ^ (f) \s~ss = for any point on the screen (Dirichlet problem). The Neumann problem can 
be treated in a similar way. For such a problem, we can in principle always divide the field into an incident one, called 
^inc {t^) and existing independently of any screen, and into a scattered field (r), produced by sources in the screen. 
The problem cannot be solved without postulating some properties of the sources. A way to do this is to introduce a 
source term J (f) in the second member of the Helmoltz equation such that this term goes to zero rapidly outside of 
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the pancake volume occupied by the screen. Then, we have [V^ + /c^]4' (r) = —J{r). Imposing Sommerfeld's radiation 
condition at infinity gives us the solution 

{r)^ I J (P) G (r, r'\ Sr'. (Al) 

J pancake 

We deduce the important fact that this potentiel ^t' (r) must be an even function of z. This is consistent with the 
Kirchhoff formula applied on the surface of Fig. 1(B). Imposing the condition 5"' [x, y, z) — 4"' {x, y, —z) implies 



^{r) = - G [f, r'j z ■ V'*' [r'j dS' (A2) 

which defines the source term Js {x, y) (surface density) by Js {x, y) = — \m\^^^)+ z ■ ^5*' [x, y, z). It is worth noting 
that the even character of and the field continuity in the aperture impose z ■ (x, y,z — Q) in the opening. 
In order to complete the problem, we must define the reflected field 5''' (r) produced by the sources when the plane 
screen contains no aperture. Since for z > there is no field, we must choose '^^ (x, y, z) = — il/* (a;, y, z) in this half 
plane. The requirement that the source field is an even function of z imposes Vl/'' (x, y, z) — —'5* (x, y, ~z) for z < 0. 
In this form, the problem is similar to the one described by Bouwkamp [lol | and it can be solved. The rest of the 
reasoning is similar to the one given for the Smythe formula. Identifying the Kirchhoff integral on the two different 
surfaces represented in Figs. 2(A) and 2(B), we obtain 

(r) = 2 y (P) z ■ V'G [r, r'^ dS' + (^j ^ - J ^ (p^ n' ■ V'G (r, - G (r, f'^ n' ■ V^' (r')]dS".(A3) 

As for the Smythe formula, we can use the symmetry properties of the field as well as its asymptotic behavior at 
infinity to transform Eq. (jA3|) into 



(f) = 2 ^ (r'^ z ■ V'G (f, r'^ dS' (A4) 
which is equivalent to the Rayleigh-Sommerfeld result given by Eq. ([2|). 
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Fig. 1. The problem of diffraction in electromagnetism. The incoming wave comes from the z < half-space and 
is diffracted by the aperture SS located in the plane screen S a.t z = 0. The unit vector n' = z used in the text is 
represented. 

Fig. 2. The two surfaces of integration for the application of the vectorial kirchhoff theorem. 
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